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Abstract — Q In this paper, a cognitive relay channel is consid- 
ered, and amplify-and-forward (AF) relay beamforming designs 
in the presence of an eavesdropper and a primary user are 
studied. Our objective is to optimize the performance of the 
cognitive relay beamforming system while limiting the interfer- 
ence in the direction of the primary receiver and keeping the 
transmitted signal secret from the eavesdropper. We show that 
under both total and individual power constraints, the problem 
becomes a quasiconvex optimization problem which can be solved 
by interior point methods. We also propose two sub-optimal 
null space beamforming schemes which are obtained in a more 
computationally efficient way. 

Index Terms: Amplify-and-forward relaying, cognitive radio, 
physical-layer security, relay beamforming. 

I. INTRODUCTION 

The need for the efficient use of the scarce spectrum in 
wireless applications has led to significant interest in the 
analysis of cognitive radio systems. One possible scheme 
for the operation of the cognitive radio network is to allow 
the secondary users to transmit concurrently on the same 
frequency band with the primary users as long as the resulting 
interference power at the primary receivers is kept below the 
interference temperature limit (TJ. Note that interference to 
the primary users is caused due to the broadcast nature of 
wireless transmissions, which allows the signals to be received 
by all users within the communication range. Note further 
that this broadcast nature also makes wireless communica- 
tions vulnerable to eavesdropping. The problem of secure 
transmission in the presence of an eavesdropper was first 
studied from an information-theoretic perspective in J2) where 
Wyner considered a wiretap channel model. In 0, the secrecy 
capacity is defined as the maximum achievable rate from the 
transmitter to the legitimate receiver, which can be attained 
while keeping the eavesdropper completely ignorant of the 
transmitted messages. Later, Wyner's result was extended 
to the Gaussian channel in 0. Recently, motivated by the 
importance of security in wireless applications, information- 
theoretic security has been investigated in fading multi-antenna 
and multiuser channels. For instance, cooperative relaying 
under secrecy constraints was studied in Il9l- |TPD . In ifTTIl . for 
amplify and forwad relaying scheme, not having analytical 
solutions for the optimal beamforming design under both 
total and individual power constraints, an iterative algorithm 

1 This work was supported by the National Science Foundation under Grants 
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is proposed to numerically obtain the optimal beamforming 
structure and maximize the secrecy rates. 

Although cognitive radio networks are also susceptible to 
eavesdropping, the combination of cognitive radio channels 
and information-theoretic security has received little attention. 
Very recently, Pei et al. in lfl2ll studied secure communi- 
cation over multiple input, single output (MISO) cognitive 
radio channels. In this work, finding the secrecy-capacity- 
achieving transmit covariance matrix under joint transmit and 
interference power constraints is formulated as a quasiconvex 
optimization problem. 

In this paper, we investigate the collaborative relay beam- 
forming under secrecy constraints in the cognitive radio net- 
work. We first characterize the secrecy rate of the amplify- 
and-forward (AF) cognitive relay channel. Then, we formulate 
the beamforming optimization as a quasiconvex optimization 
problem which can be solved through convex semidefinite pro- 
gramming (SDP). Furthermore, we propose two sub-optimal 
null space beamforming schemes to reduce the computational 
complexity. 

II. Channel Model 

We consider a cognitive relay channel with a secondary user 
source S, a primary user P, a secondary user destination D, 
an eavesdropper E, and M relays {i? m }^f =1 , as depicted in 
Figure [T] We assume that there is no direct link between S 
and D, S and P, and S and E. We also assume that relays 
work synchronously to perform beamforming by multiplying 
the signals to be transmitted with complex weights {w m }. We 
denote the channel fading coefficient between S and R m by 
g m S C, the fading coefficient between R m and D by h m G C, 
R m and P by k m £ C and the fading coefficient between 
R m and E by z m S C. In this model, the source S tries 
to transmit confidential messages to D with the help of the 
relays on the same band as the primary user's while keeping 
the interference on the primary user below some predefined 
interference temperature limit and keeping the eavesdropper E 
ignorant of the information. It's obvious that our channel is a 
two-hop relay network. In the first hop, the source S transmits 
x s to relays with power i?[|x s | 2 ] = P s . The received signal at 
the m th relay R m is given by 



Vr.ra — Qm^s i f]n 



(1) 




Fig. 1. Channel Model 



where rj m is the background noise that has a Gaussian distri- 
bution with zero mean and variance of N m . 

In the AF scenario, the received signal at R m is directly 
multiplied by l m w m without decoding, and forwarded to D. 
The relay output can be written as 



The scaling factor, 



^\g m \ 2 P s + N„ 



(2) 



(3) 



. There are two kinds of 



is used to ensure ^[lav,™! 2 ] = | 
power constraints for relays. First one is a total relay power 
constraint in the following form: ||w|| 2 = w'w < Pt where 
w = [wi,...wm] T and Pt is the maximum total power. 
(■) T and (-)t denote the transpose and conjugate transpose, 
respectively, of a matrix or vector. In a multiuser network such 
as the relay system we study in this paper, it is practically more 
relevant to consider individual power constraints as wireless 
nodes generally operate under such limitations. Motivated by 



this, we can impose 



< p m Vm or equivalently |w| < p 



where | • | 2 denotes the element-wise norm-square operation 
and p is a column vector that contains the components {p m }- 
p m is the maximum power for the m th relay node. 

The received signals at the destination D and eavesdropper 
E are the superposition of the messages sent by the relays. 
These received signals are expressed, respectively, as 

M 

Vd = ^ h m w m l m 

(gm%s + 77 m ) + no, and (4) 

m— 1 
M 



m— 1 



(5) 



where no and m are the Gaussian background noise com- 
ponents with zero mean and variance iVo, at D and E, 



respectively. It is easy to compute the received SNR at D 
and E as 



2~^m= 1 ^ J vaQm}'ra^m\ Ps 



V M \h 1272 \ w I2JU 



and 
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The secrecy rate is now given by 

R s = I(x s ;y d ) - I(x s ;y e ) 
= log(l + T d ) - log(l + r e ) 



(6) 
(7) 



(8) 
(9) 



log 



V \z l 2 Z 2 

Z^m=l \ 4 "m | l„ 



N< 



o 



V M I/, 1272 ly, 12 7V 



A/ 



A 

2/2 



m =i i^ml 2 ^^™! 2 ^™ + iv 



Em=l 
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(10) 



where /(•; ■) denotes the mutual information. The interference 
at the primary user is 



M 



M 



\w m \ N m . 



(11) 
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In this paper, under the assumption that the relays have 
perfect channel side information (CSI), we address the joint 
optimization of {w m } and hence identify the optimum col- 
laborative relay beamforming (CRB) direction that maximizes 
the secrecy rate in (fTob while maintaining the interference on 
the primary user under a certain threshold, i.e,. A < 7, where 
7 is the interference temperature limit. 

III. Optimal Beamforming 



Let us define 






h g = [K9ih, 


••> ^*m9m^m] t 'j 


(12) 


h z = [z*9ih, ■ 
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(13) 


h k = [k*9ih, 


■•i^mSm'm] , 


(14) 


D h = EMag(|fci 


2 Z 2 A 1 ,...,|/ 1m | 2 4^m), 


(15) 


D z = Diagflzi 


2 IIN 1 ,...,\zm\ 2 1 2 m N m ), and 


(16) 


D k = Diagflfci 
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(17) 



where superscript * denotes conjugate operation. Then, the 
received SNR at the destination and eavesdropper, and the 
interference on primary user can be written, respectively, as 



= Pswthgh^w 
w^DhW + Nq ' 

_ P s w t h z h z T w 
e ~~ wtD z w + N ' 
A = Pswthkhk'w + w^Dtw. 



(18) 

(19) 
(20) 



With these notations, we can write the objective function of 
the optimization problem (i.e., the term inside the logarithm 
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i (O) as 

i + r rf i 
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wtD z w+W 

w'DhW + N a + P s w t h g hg 1 w w^T> z w + N 

wtD z w + N + P s xv*ih. z hjw wtD h w + Af 

N + *r((D h + P«hghg + )wwt) JV + tr(D z wwt) 



(21) 



N + tr((D z + P s h z h z f )wwt) N + tr(D h wwt) ' 

If we denote t x = ^f^±^^ll , t 2 = 

^°~?~w^ zWW -u . define X = ww', and employ the semidef- 
inite relaxation approach, we can express the beamforming 
optimization problem as 

max t\t2 

X,ti,t 2 

s.t tr (X^h + Pshghg 1 -fi (D z + P s h z h z t ))) >JV (ti 
tr (X (D z - t 2 D h )) > JV (t 2 -l) 
tr (xfDk + Pshkhk^) < 7 



A. Beamforming in the Null Space of Eavesdropper's Channel 
(BNE) 

We choose w to lie in the null space of h z . With this as- 
sumption, we eliminate P's capability of eavesdropping on D. 
Mathematically, this is equivalent to | 5Z m=1 z m g m l m w m \ 2 = 
h z ^w| 2 = 0, which means w is in the null space of h z ^. We 
can write w = H^v, where denotes the projection matrix 
onto the null space of h z ^. Specifically, the columns of Hjr 
are orthonormal vectors which form the basis of the null space 
of h.J. In our case, is an M x (M — 1) matrix. The total 
power constraint becomes w'w = v^H^H^v = < Pp. 
The individual power constraint becomes |H^v| 2 < p 

Under the above null space beamforming assumption, T e is 
zero. Hence, we only need to maximize Yd to get the highest 
achievable secrecy rate. Td is now expressed as 



and diag(X.) < p, (and/or tr(X) < Py) and X y 0. 

(22) 

The optimization problem here is similar to that in JTTJ. The 
only difference is that we have an additional constraint due 
to the interference limitation. Thus, we can use the same 
optimization framework. The optimal beamforming solution 
that maximizes the secrecy rate in the cognitive relay channel 
can be obtained by using semidefinite programming with a 
two dimensional search for both total and individual power 
constraints. For simulation, one can use the well-developed 
interior point method based package SeDuMi lfT4l . which 
produces a feasibility certificate if the problem is feasible, 
and its popular interface Yalmip [15]. It is important to note 
that we should have the optimal X to be of rank-one to 
determine the beamforming vector. While proving analytically 
the existence of a rank-one solution for the above optimization 
problem seems to be a difficult taslQ, we would like to 
emphasize that the solutions are rank-one in our simulations. 
Thus, our numerical result are tight. Also, even in the case we 
encounter a solution with rank higher than one, the Gaussian 
randomization technique is practically proven to be effective 
in finding a feasible, rank-one approximate solution of the 
original problem. Details can be found in (8). 

IV. Sub-Optimal Null Space Beamforming 

Obtaining the optimal solution requires significant compu- 
tation. To simplify the analysis, we propose suboptimal null 
space beamforming techniques in this section . 

2 Since we in general have more than two linear constraints depending on the 
number of relay nodes and since we cannot assume that we have channels 
with real and positive coefficients, the techniques that are used in several 
studies to prove the existence of a rank-one solution (see e.g., [5|, |8 |,and 
references therein) are not directly applicable to our setting. 
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vtH^ T D h Hf v + A^o 
The interference on the primary user can be written as 

A = P^vtH^hkhk^v + v'H^DkH^v. (24) 

Defining X = vv, we can express the optimization problem 
as 

max t 
x,t 

s.t tr (X (PsH^hghg 1 !!^ - tH^B h U^ > N t 

tr (X (H^D k H^ + P s lij ^khk+H^) < 7 

and diag(H^XH^) < p, (and/or tr(X) < P T ) and X t 0. 

(25) 

This problem can be easily solved by semidefinite program- 
ming with bisection search iflOl . 

B. Beamforming in the Null Space of Eavesdropper's and 
Primary User's Channels (BNEP) 

In this section, we choose w to lie in the null space 
of h z and hk- Mathematically, this is equivalent to re- 
quiring \Y,m=i z m9mlmW m \ 2 = |h z f w| 2 = 0, and 
\Y,m=i k rn9ml m w m \ 2 = |hk f w| 2 = 0. We can write w = 



H^ fe v, where Hj; k . denotes the projection matrix onto the 
null space of h.J and hk^. Specifically, the columns of H^ fc 
are orthonormal vectors which form the basis of the null space. 
In our case, H^ fe is an M x (M — 2) matrix. The total power 

constraint becomes w^w = vtH^ fe ^H^ fc v = v^v < Py. 
The individual power constraint becomes |H^- fe v| 2 < p. 

With this beamforming strategy, we again have T e = 0. 
Moreover, the interference on the primary user is now reduced 
to 



A : 



E 

rn—l 



\k m \ 2 lL\w m \ 2 Nr, 



vtH^DkHkv (26) 
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which is the sum of the forwarded additive noise components 
present at the relays. Now, the optimization problem becomes 



max t 
x,t 



s.t tr fx fp.HiJhghtHifc - miJr> h ni k )) > n ( 



<r(x(H^ f D k H^)) < 7 

and diag(H^XH^ f ) < p, (and/or tr(X) < P T ) 
and X >- 0. 



(27) 



Again, this problem can be solved through semidefinite pro- 
gramming. With the following assumptions, we can also obtain 
a closed-form characterization of the beamforming structure. 
Since the interference experienced by the primary user consists 
of the forwarded noise components, we can assume that the 
interference constraint A < 7 is inactive unless 7 is very small. 
With this assumption, we can drop this constraint. If we further 
assume that the relays operate under the total power constraint 
expressed as v'v < Pp, we can get the following closed-form 
solution: 



max 

vtv<P t 



max 



P^vtH^Jh^H-L, 



vtv<p t vtH- 1 , f D h HJ 



An 



v+v<P t v f 



— PA, 



Pt 1 



Aoj 

Pt 



where A max (A, B) is the largest generalized eigenvalue 
of the matrix pair (A, B) 0. Hence, the maximum se- 
crecy rate is achieved by the beamforming vector v opt = 
<;u where u is the eigenvector that corresponds to 
A max (H^hgh^H^H^DhH^ + and <r is chosen 

to ensure v f opt v opt = P T . 

V. Multiple Primary Users and Eavesdroppers 

The discussion in Section [Hi] can be easily extended to 
the case of more than one primary user in the network. 
Each primary user will introduce an interference constraint 
Ti < 7i which can be straightforwardly included into d22l i. 
The beamforming optimization is still a semidefinite pro- 
gramming problem. On the other hand, the results in Section 
IIITl cannot be easily extended to the multiple-eavesdropper 
scenario. In this case, the secrecy rate for AF relaying is 
R s = I(x s ;yd) — maxi I(x s ; y e ,i), where the maximization 
is over the rates achieved over the links between the relays 
and different eavesdroppers. Hence, we have to consider the 
eavesdropper with the strongest channel. In this scenario, the 
objective function cannot be expressed in the form given in 

3 For a Hermitian matrix A 6 c«x« anc j positive definite matrix B £ 
n , (A, i/>) is referred to as a generalized eigenvalue - eigenvector pair of 



(A, B) if (A, ip) satisfy Atp = XBtp fl3l . 
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Fig. 2. AF secrecy rate vs 
10,7 = 0dS. 
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([Tol l and the optimization framework provided in Section [TTT] 
does not directly apply to the multi-eavesdropper model. 

However, the null space beamforming schemes discussed in 
Section [IV] can be extended to the case of multiple primary 
users and eavesdroppers under the condition that the number of 
relay nodes is greater than the number of eavesdroppers or the 
total number of eavesdroppers and primary users depending 
on which null space beamforming is used. The reason for this 
condition is to make sure the projection matrix exists. 
Note that the null space of i channels in general has the 
dimension M x (M — i) where M is the number of relays. 

VI. Numerical Results and Discussion 

We assume that {g m }, {h m }, {z m }, {k m } are complex, 
circularly symmetric Gaussian random variables with zero 
mean and variances a 2 g , <r\, a\ and a\ respectively. In this 
section, each figure is plotted for fixed realizations of the 
Gaussian channel coefficients. Hence, the secrecy rates in the 
plots are instantaneous secrecy rates. 

In Fig. |2 we plot the optimal secrecy rates for the amplify- 
and-forward collaborative relay beamforming system under 
both individual and total power constraints. We also pro- 
vide, for comparison, the secrecy rates attained by using 
the suboptimal beamforming schemes. The fixed parameters 
are a g = 10, cr^ = l,cr z = 1,07s = 1, 7 = OdB, and 
M = 10. Since AF secrecy rates depend on both the source 
and relay powers, the rate curves are plotted as a function 
of Pp/P- We assume that the relays have equal powers in 
the case in which individual power constraints are imposed, 
i.e., pi = Pt/M. It is immediately seen from the figure 
that the suboptimal null space beamforming achievable rates 
under both total and individual power constraints are very 
close to the corresponding optimal ones. Especially, they are 
nearly identical in the high SNR regime, which suggests that 
null space beamforming is optimal at high SNRs. Thus, null 
space beamforming schemes are good alternatives as they 
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Fig. 3. AF secrecy rate vs. Pt/Ps- a g = 10, ah = l,<Tz = 2, a^ = 
4, M = 10, 7 = lOcLB. 
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Fig. 4. AF secrecy rate vs. interference temperature 7. cr g = 10, ah = 
2, a z = 2, cr fe = 4, M = 10, P s = P T = OdB. 



are obtained with much less computational burden. Moreover, 
we interestingly observe that imposing individual relay power 
constraints leads to small losses in the secrecy rates. 

In Fig. [3] we change the parameters to a g = 10, 07, = 
l,a z = 2,a k = 4, 7 = IQdB and M = 10. In this case, 
channels between the relays and the eavesdropper and between 
the relays and the primary-user are on average stronger than 
the channels between the relays and the destination. We note 
that beamforming schemes can still attain good performance 
and we observe similar trends as before. 

In Fig. |U we plot the optimal secrecy rate and the se- 
crecy rates of the two suboptimal null space beamforming 
schemes (under both total and individual power constraints) 
as a function of the interference temperature limit 7. We 
assume that Pt = P s = OdB. It is observed that the secrecy 



rate achieved by beamforming in the null space of both the 
eavesdropper's and primary user's channels (BNEP) is almost 
insensitive to different interference temperature limits when 
7 > —AdB since it always forces the signal interference to 
be zero regardless of the value of 7. It is further observed 
that beamforming in the null space of the eavesdropper's 
channel (BNE) always achieves near optimal performance 
regardless the value of 7 under both total and individual power 
constraints. 

VII. Conclusion 

In this paper, collaborative relay beamforming in cognitive 
radio networks is studied under secrecy constraints. Optimal 
beamforming designs that maximize secrecy rates are investi- 
gated under both total and individual relay power constraints. 
We have formulated the problem as a semidefinite program- 
ming problem and provided an optimization framework. In 
addition, we have proposed two sub-optimal null space beam- 
forming schemes to simplify the computation. Finally, we have 
provided numerical results to illustrate the performances of 
different beamforming schemes. 
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